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ABSTRACT 

We have studied the wisps to the northwest of the Crab pulsar as part of a multi- 
wavelength campaign in the visible and in X-rays. Optical observations were obtained 
using the Nordic Optical Telescope in La Palma and X-ray observations were made 
with the Chandra X-ray Observatory. The observing campaign took place from Oc- 
tober 2010 until September 2012. About once per year we observe wisps forming and 
peeling off from (or near) the region commonly associated with the termination shock 
of the pulsar wind. We find that the exact locations of the northwestern wisps in the 
optical and in X-rays are similar but not coincident, with X-ray wisps preferentially 
located closer to the pulsar. This suggests that the optical and X-ray wisps are not 
produced by the same particle distribution. Our measurements and their implications 
are interpreted in terms of a Doppler-boosted ring model that has its origin in MHD 
modeling. While the Doppler boosting factors inferred from the X-ray wisps are con- 
sistent with current MHD simulations of PWNe, the optical boosting factors are not, 
and typically exceed values from MHD simulations by about a factor of 4. 

Key words: ISM: supernova remnants - pulsars: individual: Crab - radiation mech- 
anisms: non-thermal - ISM: supernova remnants: individual: Crab Nebula 



1 INTRODUCTION 

The Crab Nebula is one of the most studied targets in the 
sky as it is bright and observable over a very broad spectral 
range. The Nebula is a remnant from a supernova explosion 
that was observed on earth in 1054 C.E. Located at a dis- 
tance of ~2 kpc, the system is powered by a pulsar of spin- 
down luminosity L ~ 5 x 10 38 erg s _1 and period P ~ 34 ms. 
The history and general properties of the system are nicely 



summarized in the review by Hester ( 2008 ) 



Amongst the most prominent features of the inner neb- 
ula are the oppositely-directed jets, a torus and a rich vari- 
able "wisp" structure (Scargle|1969| more prominent to the 



northwest than the southeast. In the past decade, there have 
been several detailed observations in various wavelengths 
that observed the jet, found a knot at the beginning of the jet 
close to the pulsar (about 0.7"), and reconfirmed the wisps 
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(Hes ter et al.|1995||Bietenholz et al.|2004||Hester|2008| ). Ob- 
servations with the Hubble Space telescope (HST) showed 
the the optical wisps form and dissipate over timescales of 
months (Hes ter et al.| [2002). In addition, variable (in both 
position and time) X-ray-emitting knots also are present, 
especially to the southeast of the pulsar (Weisskopf et al. 
[2012t . 



The discovery of 7-ray flaring in 2010 September (Ta- 
Abdo et al. 2011) stimulated a renewed 



vani et al. 



2011 



interest in the Crab Nebula. Here we present the results of 
two observing campaigns, one in the visible and one in X- 
rays inspired by the search for the origin of the 7-ray flaring. 
In this paper we concentrate on characterizing the wisps to 
the northwest of the pulsar. As part of the campaign, we 
also obtained data during during the 7-ray flaring activity 
of April 2011 ( |Buehler et"aT|2012) . However, we see no di- 
rect correlation between the appearance of either an optical 
or an X-ray wisp and the 7-ray flaring. 



We describe the optical ({ 2.1) and X-ray (£2.2) obser- 
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Table 1. Optical observation dates and seeing resolution 



Figure 1. Optical image of the inner portion of the Crab Nebula 
taken on 13th of April 2011 and during a 7-ray flare (Buehler 
|et al.|2012| . North is up and east is to the left. The square covers 
a field of view of 50" x 50". The long, narrow rectangle at a 
position angle of 119° outlines the region used to measure the 
radial profile. Two ellipses are fitted to the two most prominent 
wisps at 8.0 arcsec and 10.1 arcsec distance from the pulsar. The 
ellipses have been anchored by fits to the peak positions in the 
northwest. The radial profiles of this image are shown in Fig. [3] 
The aspect ratio of the ellipses is 0.6 and the center of the ellipses 
are offset from the pulsar by 0.8 /; in the line of the pulsar and 
the peak emission of the wisps. Low frequency background has 
been removed from this image for better visibility. The units are 
arbitrary. 



vations. Data reduction and derived physical parameters are 
presented as follows: wisp radial profile s (j]3.1| ), wisp veloci- 
ties (£3.2) and wisp azimuthal profiles ({ 3.3). In ^4] we make 
use of a Doppler-boosted ring model to the observations and 
discuss implications. We briefly summarize our findings in 
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2 THE OBSERVATIONS 
2.1 The Optical Observations 

Optical data was obtained using the 2.56m Nordic Optical 
Telescope (NOT) located in the Observatorio del Roque de 
los Muchachos on the Canary island La Palma. Altogether 
24 images were taken from November 2010 until September 
2012 using the Andalucia Faint Object Spectrograph and 
Camera (ALFOSC). The instrument employs a 2048 x 2048 
E2V chip with a gain of 0.327e~/ADU, a readout noise of 
4.2e _ and a scale of 0.19^/pixel giving a field of view of 6.5 x 
6.5 arcmin. Two to four exposures of 200s - 300s were made 
through the I-band filter (NOT filter #12, A c around 800nm) 
at each epoch. The effective seeing resolution depends on 
observing conditions and ranges between 0.49" and 1.1". 
Table [l] lists the dates and seeing resolution for each optical 
observation. 

The images were processed first by subtracting the aver- 
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1.10 


56176.2 


2012/09/06 


0.75 



a Modified Julian Date of start of observation. 
b UTC date. 

c Seeing FWHM in arcsec. 



age bias image and then dividing by a twilight-sky flat-field 
image. The fringe pattern was removed using a fringe correc- 
tion image created from archival NOT data taken using the 
same instrument configuration. After scaling and subtract- 
ing the fringe correction image, the fringe pattern was no 
longer visible. Individual frames were then registered using 
8 stars in the vicinity of the pulsar and then averaged. 

We checked the NOT plate scale using an Hubble Space 
Telescope (HST) legacy archive image of the Crab pulsar 
obtained with the Wide Field and Planetary Camera 2 
(WFPC2) on Aug 14, 1995. Using the Planetary Camera 
(PC) portion of the image we selected two pairs of unsat- 
urated stars within 25 arcsec of the pulsar and determined 
the distance between the two stars in each pair using the 
0.05 arcsec/pix plate scale of the image. Measuring the dis- 
tances of the same pairs in the NOT image we arrived at 
NOT plate scale of 0.1902 and 0.1904 arcsec/pix. Hence the 
adopted plate scale 0.19 arcsec/pix is expected to be accu- 
rate to within 0.2%. 

Sky background was subtracted from the average frames 
by fitting two models to the "pure sky" regions around Crab: 
a constant and a low-order polynomial. The two methods 
generally agreed to better than 1.2% when extrapolated to 
the location of the pulsar. We use the constant in what fol- 
lows. Flux calibration used the comparison star sequence 
Sandberg Sollerman| ( |2009| . We used star 3 with 



given m 

I-band magnitude I = 16.26 =b 0.02 from their calibration 
sequence to determine the zero point, ZP, by fitting a Point 
Spread Function (PSF) model created from star 1 of |Sand- 



berg &; Sollerman (2009) to star 3. The conversion from sky- 
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Figure 2. This Chandra X-ray image of the inner portion of the 
Crab Nebula is the sum of the observations from the 14th to the 
28th of April 2011 during the 7-ray flare. The axes are ACIS 
detector pixels (0.492^/pixel. The best fitting ellipse of the wisp 
emission has an aspect ratio of 0.49 and an offset of 0.9" at a 
position angle of 110° (|Weisskopf et al. 2012). The color scale is 
in counts. The black feature with zero counts is the burned out 
image of the pulsar. The color scale are photon counts per ACIS 
pixel in the 0.5-8 keV band. 



subtracted counts, N (ADU), to physical fluxes, F (/i Jy per 
pixel), is then 



F = 2550* 10-°- 4 * (zp - 2 - 51ogiv) 



(1) 



Since the calibration was obtained from a star in the Crab 
image, variations in atmospheric transparency are effectively 
accounted for and removed. To check the accuracy of the cal- 
ibration we measured the total flux of a 5.9" x 7.8" wisp- free, 
rectangular region ~ 40" east of the Pulsar and of apparent 
uniform surface brightness. The rms scatter of the flux from 
this region is 3.3% over the 24 observations. We ascribe most 
of this scatter to flat-fielding-, calibration- and background- 
determination errors as the expected scatter, considering 
only the photon noise and the readout noise, is < 0.1%. 
Fig. [l] shows one of the 24 optical images. 



2.2 The X-ray Observations 

Subsequent to the discovery of 7-ray flaring in 2010 Septem- 
ber (Tavani et al. 2011 Abdo et al.|2011 ) a series of Chandra 
X-ray Observatory observations was initiated on approxi- 
mately a monthly cadence. Five Target of Opportunity ob- 
servations were performed during the April 2011 7-ray flare 
and the average of those 5 images is shown in Fig. [2] Chan- 
dra observations have continued, now as part of the Chandra 
general observer program, yielding a sequence which covers 
approximately two years. As with our optical observations, 
the coverage is not continuous, since Chandra is not allowed 
to point within 45° of the Sun. The solar restriction results 
in about a 90 day interval during the summer when the Crab 
may not be observed. A list the Chandra observations used 
for this paper is given in Table [2] 



All data were taken with the back-illuminated ACIS-S^H 
detector at the focus of the telescope in the energy band 0.5- 
8 keV. Owing to the Crab's high flux, the observations were 
taken with the shortest possible frame time of 0.2s to mini- 
mize pileup. Even with this frame time, pileup impacts the 
brighter regions of the nebula. Telemetry saturation is also 
a issue with such a bright source, so we restricted data to a 
300 x 300 ACIS-pixel (« 150" x 150") region that covers the 
bright part of the nebula. The integration time was 5 ksec 
for monitoring observations and 10 ksec for the 5 observa- 
tions following the 2011 April 7-ray flare. Due to telemetry 
saturation the effective exposure time was 600 or 1200s (Ta- 
ble^. 

For the first three observations (Table [2| the dither am- 
plitude of the spacecraft was set to zero. Since pileup charac- 
teristics may be different if the pulsar image is at the middle, 
edge or corner of a pixel, we then used a small (1" ampli- 
tude) dither for the remaining observations. This dither en- 
sures the pulsar and other sharp features in the nebula are 
averaged over a few pixels. 

Level 2 event files were created using the CIAO tool 
acis_process_events using the EDSER option ( Li et ah] 
2004) to improve the subpixel positioning. Due to severe 
pileup the pulsar does not appear in the image (Fig. |2|. 
We did not simply rely on the Chandra aspect solution, but 
re-registered images making use of the readout trail (out-of- 
time image of the pulsar) to constrain one dimension and 
the burned out image of the pulsar to constrain the other 
coordinate. 

Due to the different roll orientations of the Observatory, 
the readout trail could introduce some bias in the azimuthal 
distribution. We removed this trail by computing the aver- 
age number of counts per pixel due to the trail. Then, based 
on the actual number of counts seen in a pixel and our esti- 
mate of the number due to the readout streak, we randomly 
rejected events. Since this method is statistical, it is not per- 
fect, but it will remove the principal brightness peak in the 
data due to the trail itself. 



3 DATA ANALYSIS 

Images were analyzed to characterize the shape of the wisps 
in both the radial direction from the pulsar and angular 
extent about the pulsar. Our analysis to determine the ra- 
dial distribution through the wisp restricted data to a slim 
elongated rectangle 3" in width at a position angle of 119° 
measured positive east of north. The rectangular region is 
shown in Figs, [I] and [2] Optical and X-ray data were binned 
in radial bins of 0.11" and 0.492" respectively. The anal- 
ysis of azimuthal distributions utilized 0.5"- wide elliptical 
annuli, on the sky and centered on the peak of the radial 
distribution. 



3.1 Analysis of the radial profile of the wisps 

An example of an optical and an X-ray radial profile mea- 
sured close in time is shown Fig. [3] We see two prominent 
peaks in the optical and the presence of at least three peaks 



http://asc.harvard.edu/proposer/POG/ Chapter 6. 
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Table 2. Time-ordered list of Chandra observations with expo- 
sure time used for this paper. 



ObsID a 


MJD 6 


Date c 


Time d 


13139 


55467 2 


2010/09/28 


600 


13146 


55497.7 


2010/10/28 


600 


13147 


55528.5 


2010/11/28 


600 


13204 


55576.0 


2011/01/15 


600 


13205 


55608.7 


2011/02/16 


600 


13206 


55635.2 


2011/03/15 


600 


13207 


55663.6 


2011/04/12 


600 


13150 


55665.0 


2011/04/14 


1200 


13151 


55665.6 


2011/04/14 


1200 


13152 


55667.4 


2011/04/16 


1200 


13153 


55673.0 


2011/04/22 


1200 


13154 


55679.3 


2011/04/28 


1200 


13208 


55788.4 


2011/08/15 


600 


13209 


55819.2 


2011/09/15 


600 


13210 


55849.2 


2011/10/15 


600 


13750 


55892.7 


2011/11/27 


600 


13751 


55938.1 


2012/01/12 


600 


13752 


55967.4 


2012/02/10 


600 


14416 


56005.3 


2012/03/19 


600 


13754 


56039.1 


2012/04/22 


600 


13755 


56155.6 


2012/08/16 


600 


13756 


56181.9 


2012/09/11 


600 



a Chandra observation identifier. 
b Start of observation MJD. 
c UTC date. 

d Approximate effective exposure time in seconds. 



Table 3. Velocities of the wisps identified in Figs. [5] and [6] 
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Figure 3. The upper panel shows the radial distribution for the 
optical data for April 13, 2011 (in units of iiJy / arcsec 2 ) and the 
bottom panel the X-ray radial distribution for the average from 
April 14 until April 28, 2011 (photon counts per ACIS pixel). The 
middle curve is the same as the uppermost curve but with low 
frequency terms subtracted. The solid lines are a model fit based 
on a number of Gaussians. The vertical at 10 /; shows that the 
optical and X-ray peaks do not exactly coincide. 



that the optical wisps are brightest at this distance and then 
fade at larger distances. Fig. ^] bottom shows the same plot 
in X-rays averaged with the same sampling interval as for 
the optical. The innermost X-ray arc, at ~ 7.5", is brightest. 
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0.014 
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0.24 


box solid 


0.017 


0.19 


0.28 
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0.020 


0.23 


0.31 


star solid 
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triangle up 
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0.36 


triangle down 


0.011 


0.12 


0.19 


cross 


0.009 


0.10 


0.16 



a Apparent velocity on the sky. 
b Deprojected velocity. 

in the X-ray profile. In general, the peak closest to the pulsar 
is most usually identified with the location of the termina- 
tion shock. Fig. [5] demonstrates that the X-ray peak at 10" 
is located at a slightly different (and smaller) distance than 
the peak in the optical. 

Fig. ^]top is the radial profile of the average of all the 
optical observations. Notice the strong peak at 8" showing 



3.2 The radial evolution of the wisps as a function 
of time 

Fig. [5] shows most of the optical and X-ray radial pro- 
jections as a function of the angular distance from the 
pulsar along the position angle of 119° as a function of 
time. Low frequency terms (background) have been sub- 
tracted for better visibility. The algorithm applied was the 
TSpectrum::Background()-method from the ROOT-package 
dMorhac et al.|1997|). In the fi gure we also trace the outward 
motion of a particular peak by drawing a line through the 
position of what appears to be the same peak but at different 
times. There appear to be possibly 5 distinct progressions 
of peaks in the optical data and possibly 7 in X-rays. We 
have not attempted to connect progressions across summer 
gaps, but one could do so. As the slopes of the lines indicat- 
ing the outward progression are different, so are the inferred 
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Figure 4. The upper panel shows radial profile of for the aver- 
age of all the optical observations (in units of fijy /arcsec 2 ). The 
lower panel shows the average radial profile for all the X-ray data 
(photon counts) with the exception of ObsIDs 13150 to 13154 
to avoid biasing the average due to the concentrated coverage in 
April 2011 and to maintain a similar sampling cadence as for the 
optical data. The middle panel repeats the upper panel but with 
low frequency terms subtracted. The vertical line at 8 /x shows that 
the brightest and closest optical and X-ray peaks do not exactly 
coincide. 



velocities. These are given in Table [3] where we list the mea- 
sured quantity i.e. the velocity on the sky in "/day, and, 
assuming a distance of 2 kpc, the physical velocity and the 
deprojected physical velocity. For deprojection we use an 
observation angle of ~ 57 . 

Fig. [5] also compares the motion of the peaks in the 
radial profiles as a function of time. The figure seems to 
indicate that each time an X-ray wisp appears close to the 
pulsar, so does an optical wisp that is slightly further away 
from the pulsar. 

Fig. [7] presents a different yet equally interesting, pic- 
ture and comparison of the evolution of the peak fluxes in 
the azimuthal distributions as a function of time. To con- 
struct this figure we used a time spacing of 10 days. If an 
observation took place at any time within the 10 day interval 
the data were included as a column in the figure. For short 
gaps in the time sequence, we performed a linear average of 
the closest columns that contain data. Thus, if there was a 
two column gap, we would add two thirds of the previous 
observation to one third of the following observation to fill 
in the first missing column. For the second missing column, 
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Figure 5. Twenty-three optical and 17 X-ray radial profiles with 
low-frequency terms removed are plotted as a function of distance 
northwest of the pulsar. For clarity when one or more observa- 
tions occur very close in time, we omit them from the figure. The 
omitted observations are optical - 56140.2 and X-ray - 55665.0, 
55665.6, 55667.4, 55673.0, & 56039.1 Symbols are placed on the 
distance axis to show the position of the peak in the profile. So 
doing gives the false illusion that the symbols are offset from the 
peak value. The lines, simply best guesses, are an interpretation 
of the time evolution of the position of a particular peak. 



the weights were reversed. If the wisps in two data sets over- 
lap, this method will cause the wisp to appear to move from 
the location in the first data set to the location seen in the 
second across the gap. If the features do not overlap, then 
features in the first data set will appear to fade as new fea- 
tures appear to grow. We have determined that the gap in 
the summer break is too long for our interpolation method 
to work and thus they are empty. 

All the features seen in Fig. [5] can also be seen pictori- 
ally in Fig. [7] Note that the width of neighboring peaks in 
the X-ray and optical do not appear to be correlated. In the 
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55500 55600 55700 55800 55900 56000 56100 56200 
MJD 

Figure 6. The radial positions of the wisp peaks as a function 
of time. Red and blue are used to indicate an optical or an X-ray 
wisp respectively. The symbols are the as in same as in Fig. [5] 



X-ray portion of Fig. [7] we also see the slow outward motion 
of the outer boundary where the color changes from green to 
blue. This boundary is close to 20" from the pulsar for the 
earliest observations, yet moves to near 25" by the end of 
the sequence. This outward motion of the boundary roughly 
matches the apparent outward motion of the brighter wisps. 
Finally, some X-ray peaks do not appear to have nearby 
optical companion. As a conclusion we can state that the 
positions of the optical wisps do not align with the posi- 
tions of the x-ray wisps. That would seem to imply that the 
individual evolution of optical and X-ray wisps are different. 
Fig. [7] also makes it clear that new wisps form in the inner 
region roughly once per year. 




800 
MJD-55000. 




800 
MJD-55000. 



Figure 7. The two panels compare the radial evolution of the 
optical (upper panel) and X-ray (lower panel) wisps. For clarity, 
data were interpolated between observations but not across the 
large gaps imposed by sun constraints. 



3.3 Analysis of the azimuthal profile of the wisps 

Good statistics for the optical data allow us to measure the 
azimuthal intensity profile for each observation. The wisps 
have an ellipse-like shape in the azimuthal direction around 
the pulsar, presumably because they are formed in a ring 
more or less in the equatorial plane of the pulsar. In this 
case the ellipse-like shape is simply a consequence of pro- 
jecting a circular ring onto the plane of the sky. In X-rays 
one finds that the aspect ratio of the innermost ellipse (the 
location of the presumed termination shock) is 0.49 and the 
pulsar is displaced along the axis of symmetry by ~ 0.9" 
below the plane of the ring and along the minor axis which 
is at a position angle of 110° (Weisskopf et al. 2012). It is 
more difficult to determine these parameters in the visible 
primarily because a complete ring is not observed. Neverthe- 
less, averaging all the optical images and assuming this geo- 
metrical model, we find a similar offset for the pulsar (0.8") 
along the spin axis as in X-rays, but with an aspect ratio of 
0.6. The different aspect ratio is surely a consequence of our 
conclusion that the optical and X-ray wisps are not formed 
in the same regions, with the optical wisps further from the 
pulsar. Fig. [8] shows the azimuthal distribution of the in- 
tensity for two of the optical wisps at distances of 8.0" and 
10.1" from the pulsar. The distributions shown in the figure 



x 

2 




Figure 8. The azimuthal distributions of two optical wisps ob- 
served on the 13th of April 2011 at distances of 8.0" and 10.1" 
from the pulsar. The angle is the azimuthal angle of the depro- 
jected circle. Zero degrees was set by the peak of the distribution. 
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was based on a 0.5 ;/ -wide annulus on an ellipse with the 0.6 
aspect ratio. 

The X-ray data are noisier than the optical data so we 
considered both the average of the 2011 April data and the 
average of all except for the 2011 April data. Within errors, 
the width of these two distributions were identical. The az- 
imuthal variation in an annulus 0.5" wide, 8.0" from the 
pulsar, and centered at the position angle of 119° is shown 
in Fig. [5] Comparing the X-ray angular distribution with the 
optical, it is important to note the latter has a sharper peak 
than that seen in X-rays. The FWHM of the distribution 
in optical is around 30° while the X-ray distribution has a 
FWHM of around 70°. 



4 MODEL FITTING AND THEORETICAL 
DISCUSSION 

The MHD model developed with success in the past years 
dKomissarov fe Lyubarsky|2003[|2004||Del Zanna et al.|2004| 
Buccianti ni et al.||2005| |Del Zanna et al.||2006| |Volpi et al.| 
2008 ; [Camus et al.|2009| ) explains the existence of rings and 
wisps observed in pulsar wind nebulae (PWNe) in terms of 
axisymmetric enhancements of the emissivity immediately 
downstream of the termination shock. Due to the anisotropy 
of the pulsar wind luminosity, and the related oblate shape of 
the termination shock ( |Del Zanna et al.|2004| ) , axisymmetric 
regions form inside the nebula where the emissivity toward 
the observer is higher, either because of local compression of 
the magnetic field, or because of Doppler boosting. In gen- 
eral it is found that typical flow speeds immediately down- 
stream of the termination shock can be as high as 0.7c, and 
the emission is dominated by Doppler boosting. 

Komissarov & Lyubarskyl (12004) have shown that both 



the wisps and the knot observed in the Crab Nebula can be 



qualitatively explained by the geometry of the flow (Fig. 10 ). 
In particular the association of the wisps with the wind ter- 
mination shock is due to the fact that they originate in the 
flow immediately downstream of the shock itself. A simpli- 
fied model for a wisp can be built on the assumption that a 
wisp is due to a torus- or ring-like region within the nebula. 
This region is bright or dim depending on Doppler boosting 
due to the relative direction of the observer, and the par- 
ticle flow. This boosted-ring model appears in alternative 
explanations for the origin of the wisps, e.g. in term of ion- 
cyclotron compression ( Spit kovsky &; Arons|2004"| or cooling 
instabilities (Foy & Hester 2009 ). The boosting is required 



in all of these models to explain the azimuthal luminosity 
profile of the wisps, and the fact that one side of the nebula 
(the front side) is brighter than the back side. This idea has 



been adopted in the past ( 


Romani & Ng 2003 Ng k Romani 


2004 Romani et al.||2005 


) to model the jet torus structure 



observed in X-rays in several PWNe, and to estimate typical 
flow speeds in these systems. 

With reference to Fig. |10| a wisp originates from a re- 
gion shaped like a torus, with major radius R Q and a minor 
radius r Q . The ratio r /R we call the thickness of the torus. 
The plane of the torus has an inclination a with respect to 
the line of sight, which, within the axisymmetric approxi- 
mation, is equal to the angle between the nebular axis and 
the plane of the sky. The fluid in the torus has a uniform 
flow speed Vfi, confined in meridianal planes (planes con- 
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Table 4. Fitting parameters for the optical wisps. It lists the 
MJD, distance to the pulsar, the flow velocity, the boosting four- 
velocity, and the minimum and maximum flow velocities still com- 
patible with the fit and the thickness of the torus. In all cases we 
found a = (3 = 37° ±4°. 



taining the nebular axis). The azimuthal component of this 
flow speed is assumed to be 0. Vfi forms an angle f3 with 
respect to the plane of the torus. 

We want to stress here that Vfi is the flow speed of 
particles in the wisp, and should not be confused with the 
observed wisp velocity, which we have discussed previously. 
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Figure 9. Left panel: The average azimuthal X-ray intensity distribution (normalized to its maximum) in an elliptical annulus at distance 
of 8 /; to the northwest of the pulsar. The dashed line is a simulated azimuthal distribution as discussed in ^] The dotted line is an 
assumed background level. Right panel: Same as the left panel but for the average of the optical data nearest the pulsar. 



The two are not the same simply because the wisp velocity 
traces the propagation of the ring in the nebula as a wave, 
and not the bulk motion of particles. For simplicity we as- 
sume uniform synchrotron emissivity inside the torus. In this 
picture the torus contains a fully ordered toroidal magnetic 
field, but models can also be built with a chaotic magnetic 
field (see below for a discussion) . The emissivity toward the 
observer is computed as in previous wor ks ([K omissarov & 
Lyubarsky||2003l [2004) |Del Zanna et al.||2QQ4| |Bucciantini 
et al [2005 MDel Zanna et al.|2006| |Volpi et al.1|2008| |Camus| 
et al.|[2009| and a synchrotron map is built integrating the 
contribution of the various parts of the torus along the line of 
sight. The emitting particles are assumed to be distributed 
as a power-law in energy with an index —2.35 in agreement 
with values estimated from optical observations and spec- 
tral modeling of the Crab Nebula ( |Bucciantini et alT]|2011 ). 
This also applies for the X-ray wisp(s). In the MHD model 
of PWNe, the wisps are assumed to be located immediately 
downstream of the termination shock, where particles are 
accelerated and magnetic field is compressed. Given that 
the wisps originate from the termination shock, synchrotron 
cooling should be negligible and the particle distribution un- 
affected. 

The wisps appear brighter northwest of the pulsar and 
almost fade to the background level in the southeast mostly 
due to Doppler boosting, given that the flow velocity is gen- 
erally directed toward the observer in the northwest, while 
on to the southeast it points away. We want to stress here 
that the plane of the torus does not necessarily need to co- 
incide (even if it is parallel) with the equatorial plane of the 
pulsar, but there can be an offset between the two. Such an 
offset, however, has no consequences in the model for the 
wisp emissivity. 

We fit the optical data to this model as follows: 



• We select a set of input parameters r /R , /3, and 
V f i. 

• An emission map is built (the emission is normalized to 
the maximum). 



• The emission map is convolved with the point spread 
function of the observation. 

• Emission profiles are extracted both along the wisp and 
along the axis of the nebula. 

• The simulated profiles are compared with the 
background-subtracted data, both along the axis and along 
the wisps. 

• This procedure is repeated until a set of input param- 
eters are found that provide a reasonable fit. The best fit 
parameters are obtained by minimizing the residuals in the 
brightest part of the wisp, where the background subtracted 
intensity is > 50% of its maximum value. 

The results for the optical wisps are presented in Ta- 
ble [4] We can also roughly quantify the uncertainty in the 
various parameters. The upper (lower) limit on the flow 
boosting velocity Vfi given in Tab. [4] corresponds to a model 
which gives a wisp whose FWHM is 15% narrower (larger) 
than the observed value. Fig. |ll| compares one of the optical 
observations and the result from a simulated map, using the 
best fit input parameters from Table [4] 

We stress here that our procedure does not take into 
account time-of-flight delays. However, given that the extent 
of the optical wisps are usually within ±20° from the axis, 
this effect should be negligible. 

There are several interesting points to note: 

• Typical flow speeds inferred from our model range be- 
tween Vfi = 0.8c and Vfi = 0.95c. Using a chaotic magnetic 
field distribution leads to even higher velocities. 

• The angle f3 is almost equal to the angle a. 

• The wisps appear to be unresolved and the ratio r /R 
is less than 0.1. 

The high values of the flow speed correspond to boost- 
ing factors that are in excess of typical values found in MHD 
simulations by about a factor between 2 and 5 depending on 
which value of /3 one chooses. This is is related to the nar- 
row extent of the optical wisps which fade to the background 
level within ±20° from the axis. This is also connected to 
the fact that ($ w a, implying that the flow speed on the 
front side of the nebula is directed toward the observer. If 
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Figure 10. Left panel: Schematic structure of the termination shock and wisps region. The pulsar is located in the central circle from 
which a relativistic radial wind originates. This wind slows down and is diverted in a termination shock with an oblique shape due to 
wind anisotropy. The large arrow indicates the direction toward the observer, while the dashed line represents the plane of the sky. The 
symbol a is the angle between the nebular axis and the plane of the sky and (3 is the angle between the plane of the wisps and the flow 
direction in the wisp's regions. Right panel: Schematic of the wisp modeled as an emitting ring. The dashed line is the plane of the sky 
and R and r Q are the major and minor radii of the torus respectively. 



f3 differs from a by more than 3° — 5°, it is not possible to 
reproduce the azimuthal wisp luminosity, independently of 
Vfi. The minimum value of Vfi is obtained for f3 — a. More- 
over, it appears that the optical wisps are consistent with 
narrow features, possibly close to emitting sheets. It is also 
interesting to note that there is a trend in the observed flow 
speeds as a function of wisp location, with higher values for 
inner wisps, as shown in Fig. |12| 

A comparison between the azimuthal profile in optical 
and X-rays has also been done. However, for the X-ray data, 
in order to have sufficient statistics, instead of a single epoch, 
we consider first an average over the entire set of observa- 
tions and only for the inner (brighter) X-ray wisp region 
located ~ 8" away from the pulsar (Fig. [9]). We also re- 
peated this analysis averaging just the 2011 April data with 
similar results. 

The same was done for the optical data, again taking 
the average over all of our epochs. Results were shown in 
Fig. [9] The resulting model parameters are r /R = 0.1, 
a = /3 = 35° both for optical and X-ray, while Vfi = 
0.6c ± 0.1c for the X-ray and V n = 0.91c ± 0.03c for the 
optical. We may associate this with a boosting four- velocity 
U f i = (V f i/c)/y/[l - (Vfi/c) 2 } of 0.75 for X-ray and 2.35 for 
optical (i.e. a factor of 4). Thus, we can safely conclude that 
the azimuthal extent of this X-ray wisp is larger than in the 
optical, and consequently the lower inferred speeds in X-ray 
~ 0.5 — 0.7c are definitely consistent with expectations from 
MHD. 

The fitting parameters for the optical wisps, however, 
raise the question if the boosted ring model is in fact ap- 
propriate. The MHD model has been developed mostly to 
reproduce Chandra observations. It is also evident from our 
study that X-ray and optical wisps are not produced by the 



same particle distribution: they do not coincide in location 
or in terms of the degree of Doppler boosting. 

It is not our intention here develop and provide an al- 
ternative explanation for the azimuthal luminosity of the 
optical wisps, but we can point to a few possible solutions: 

• Optical wisps appear to originate in the transition layer 
at the termination shock itself. This is consistent with the 
thin-sheet hypothesis and the idea that, at the shock itself, 
the flow speed is higher than in the downstream region. It 
remains to be proven that the geometry of the shock can 
allow for such configurations. 

• Emissivity is enhanced along the axis of the nebula be- 
cause of local dumpiness. Clumps are observed in X-rays, 
called by some sprites, and they have a more or less stable 
location. However, why in optical, these clumps should be 
located along the axis of the nebula is not clear. 

• If the particle distribution in a wisp is non-isotropic, 
then the emissivity has a stronger dependence on the incli- 
nation of the magnetic field with respect to the line of sight 
than in the isotropic case. It is not clear, however, what the 
degree of anisotropy needs to be and if it can be produced 
by known acceleration mechanisms. 



5 CONCLUSION 

We find that approximately once per year optical and X-ray 
wisps appear and peel off from the region commonly associ- 
ated with the termination shock of the pulsar wind. More- 
over, the X-ray and optical wisps are physically separated. 
Interestingly enough, the time interval for 7-ray flares is also 
about 1 year leading to the speculation that the flaring may 
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Figure 12. Boosting four- velocity (Vfi/c)/ — (Vfi/c) 2 ] for 
the optical wisps as a function of their distance from the pulsar. 
The square denotes the boosting four-velocity obtained from the 
averaged optical profile. Note that the data points at radii ^ 7" 
are most likely not due to wisps, but other optical features such 
as the halo ( |Hester|200"8| >. 
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Figure 11. Upper panel: Simulated synchrotron map of the op- 
tical wisp for MJD=55532 and at a distance of 7.7" derived using 
the fitting parameters in Tab. ^] The axes are in arbitrary units 
normalized to the wisp major radius R Q . Colors indicate the level 
of the flux, normalized to the maximum. Middle panel: the solid 
line is the optical profile of the wisp along the nebular axis, while 
the dashed line is obtained from the simulated synchrotron map. 
The dotted line is the (assumed) background. Lower panel: Same 
as the middle panel but for the azimuthal profile. 



somehow be associated with wisp formation. However, since 
the wisps did not form in coincidence with the 7-ray flare of 
2011 April which took place during this study, the specula- 
tion requires a theory that has the formation of the wisps 
substantially (^month) preceding a 7-ray flare. 

These wisps propagate outward with velocities pro- 
jected onto the sky ranging from 0.1 (v/c) to 0.4(v/c) (Ta- 
ble [3}. If anything, there appears to be a trend for the pro- 
jected velocity of a wisp increasing, the further it is from the 
pulsar, perhaps implying a re-acceleration mechanism. More 
likely, this apparent behavior may be a result of the more 
complicated three-dimensional geometry and there may be 
no need for re-acceleration when properly deprojected. An- 
other word of caution is that one might be confusing optical 
features in the counter-jet with structure in the wisps. 

Within the context of an MHD model, we also find that 
optical wisps are more strongly Doppler-boosted than the X- 
ray wisps. In particular, we found that the azimuthal lumi- 
nosity profile of the X-ray wisps is fully compatible with typ- 
ical boosting factors found in MHD simulations of PWNe. 
Instead, the azimuthal luminosity profile of the optical wisps 
requires particles velocities that are incompatible with the 



results of global numerical modeling of PWNe (Komissarov 
& Lyubarsky 2003|j Bucciantini et al.|2005 ). It is also inter- 



esting to note that the optical wisps as well as the X-ray 
wisps are Doppler boosted in the direction of the observer, 
which implies that the movement of particles has a compo- 
nent in the direction of the observer. This may be simply 
an accident or due to a selection effect as we tend to ob- 
serve high-flux regions and the Doppler boosting enhances 
the flux. This should be investigated with future modeling. 
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